UFIFT-HEP-95-6, Fermilab-Pub-95 /080-A 



The velocity peaks in the cold 
dark matter spectrum on Earth 

P. Sikivie^'^, I. I. Tkachev^''^ and Yun Wang^ 
^Physics Department, University of Florida, Gainesville, FL 32611 

Institute for Theoretical Physics, Santa Barbara, CA 93106 
^NASA/Fermilab Astrophysics Center, FNAL, Batavia, IL 60510 
^Institute for Nuclear Research, Russian Academy of Sciences, Moscow 117312, Russia 

(April 13, 1995) 

Abstract 

The cold dark matter spectrum on earth is expected to have peaks in velocity 

space. We obtain estimates for the sizes and locations of these peaks. To this 

end we have generalized the secondary infall model of galactic halo formation 

to include angular momentum of the dark matter particles. This new model 

is still spherically symmetric and it has self-similar solutions. Our results are 

relevant to direct dark matter search experiments. 
PACS numbers: 98.80.Cq, 95.35.+d, 14.80.Mz, 14.80.Ly 
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There is considerable evidence that galaxies are surrounded by dark halos which con- 
tribute 90% or more of the galactic mass [jl|. Identifying the nature of the dark matter is the 
goal of considerable theoretical and experimental effort. The leading dark matter candidates 
are baryons, neutrinos, weakly interacting massive particles (WIMPs) and axions. From the 
point of view of galaxy formation, axions and WIMPs have identical properties which earn 
them the name of cold dark matter (CDM). Unlike baryonic dark matter, CDM is garanteed 
to have negligible interactions other than through its gravitational effects. Unlike neutrinos, 
CDM has negligibly small primordial velocity dispersion. Studies of large scale structure 
formation support the view that the dominant fraction of dark matter is CDM. Moreover, 
if some fraction of the dark matter is CDM, it necessarily contributes to galactic halos by 
falling into the gravitational wells of galaxies. The main focus of our paper is to study the 
growth of galactic halos by the infall of cold dark matter, including the effect of angular 
momentum, and thence to derive properties of the spectrum of dark matter particles on 
Earth. 

We are motivated in part by the prospect that a direct search experiment may some 
day measure that spectrum. In particular, if a signal is found in the cavity detector of 
galactic halo axions it will be possible to measure the spectrum with great precision 
and resolution. Naturally, the question arises what can be learned about our galaxy from 
analyzing such a signal. Moreover, knowledge of the spectrum, which is identical for both 
CDM candidates, may help in the discovery of a signal. 

In many past discussions of dark matter detection on earth, it has been assumed that the 
dark matter particles have an isothermal distribution. Thermalization has been argued to 
be the result of a period of "violent relaxation" |^ following the collapse of the protogalaxy. 
If it is strictly true that the velocity distribution of the dark matter particles is isothermal, 
which seems to be a very strong assumption, then the only information that can be gained 
from its observation is the corresponding virial velocity and our own velocity with respect 
to its standard of rest. If, on the other hand, the thermalization is incomplete, a signal in a 
dark matter detector may yield additional information. 
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J. R. Ipser and one of us discussed 0] the extent to which cold dark particles are thermal- 
ized in a galactic halo and concluded that there are substantial deviations from a thermal 
distribution in that the highest energy particles have narrow peaks in velocity space. There 
is one velocity peak on earth due to particles falling onto the galaxy for the first time, one 
peak due to particles falling out of the galaxy for the first time, one peak due to particles 
falling into the galaxy for the second time, and so on. A simple topological argument shows 
that these peaks exist no matter what is the distribution of angular momenta the dark 
matter particles have with respect to the galactic center, although angular momentum may 
reduce the sizes of the peaks. The peaks due to particles which have fallen in and out of the 
galaxy a large number of times in the past are washed out by scattering in the gravitational 
wells of stars, globular clusters and large molecular clouds. Those particles are in effect 
thermalized. But the peaks due to particles which have fallen in and out of the galaxy only 
a few times in the past are not washed out by scattering. 

The width of the peaks due to the primordial velocity dispersion of the CDM particles 
is very small: AE/m ~ 10^^^ for axions and !S.Ejm ~ 10^^^ for WIMPs, where m is the 
mass of the particles and !S.E their energy spread. As just mentioned, the peaks fragment 
because of the gravitational scattering of the particles by inhomogeneities in the galaxy. 
They also fragment because of structure formation on scales smaller than that of the galaxy 
as a whole. We expect that such structure formation in the CDM component of the matter 
density would be inhibited by the tidal forces of the galaxy's gravitational field, but some of 
it no doubt occurs. At any rate, it is clear that the peaks due to particles which have fallen 
in and out of the galaxy only a small number of times in the past are not entirely washed out 
by either scattering or small structure formation. It is the purpose of the present paper to 
estimate the sizes and locations of these peaks and we will no longer concern ourselves with 
their widths. Let us point out however that the sensitivity of the search for galactic halo 
axions using the cavity detector may be increased by looking for narrow peaks. Specifically, 
it is found that the signal to noise ratio of the upcoming axion search at LLNL [Q] is increased 
by a factor I8O/1 by looking for narrow peaks, where /i is the fraction of the local axion 
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density in the largest of the peaks of width Ai? < 10 m. 

The tool we use to obtain estimates of the locations and sizes of the highest energy peaks 



makes the following assumptions: 

1) the dark matter is non-dissipative; 2) it has negligible initial velocity dispersion; 3) 
the gravitational potential of the galaxy is spherically symmetric and is dominated by the 
dark matter contribution; 4) the dark matter particles have zero angular momentum and 
therefore move on radial orbits through the galactic center. 

Below, we will generalize the model to rid it of the fourth assumption. For clarity, we 
refer to the model with the fourth assumption included as the radial infall model. 

An initial overdensity profile SMi{r) is assumed. The equations of motion for the radial 
coordinate r{a,t) of each spherical shell {a is a shell label, t is time) in the gravitational 
potential due to all the other shells must then be solved for initial conditions given by 
the Hubble expansion at some arbitrarily chosen but early time tf. r{a,ti) = H{ti)r{a,ti). 
Much progress in the analysis of the model came about as a result of the realization that 
the evolution of the galactic halo is self-similar ^ provided the initial overdensity has the 
following scale-free form: 



where Mi and 6Mi are respectively the mass and excess mass interior to rj at the initial time 
U, i.e. 



where p{ti) = 3H{ti)^ /SnG = l/GnGt"^ is the average density in a critical {Q = 1) uni- 
verse, ti having been chosen during the matter-dominated epoch. Mq and e are constants 
characterizing the model. Self-similarity means that the phase-distribution of dark matter 
particles is time-independent after all distances have been rescaled by the overall size R{t) of 
the galactic halo, all masses by the mass M{t) interior to the radius R{t), and all velocities 



is the secondary infall model of galactic halo formation 0. In its original form, this model 




(1) 



Miin) = —piuyf + SM^in) 



(2) 
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by G M (t) / R{t) . For the sake of definiteness, R{t) is taken to be the "turn-around" radius 
at time t, i.e. the radius at which particles have zero radial velocity for the first time in 
their history (see Fig.l). In a self-similar solution, the motion of each shell a is the same 
after the appropriate rescaling, i.e. r{a,t) = r,,(a;)A(t/t*(a;)), where r^{a) and t*(a) are 
respectively the turn-around radius and time of shell a. Also the mass-profile of the halo is 
time-independent after rescaling, i.e. M{r,t) = M{t)M.{r / R{t)). A(r) and M.{i) are func- 
tions of a single variable which can be accurately obtained by numerical integration. In the 
limit ^ ^ it was shown analytically that, oc ^ if < e < 2/3 and oc ^^/{i+ac) 

if 2/3 < e < 1. The numerical integrations confirm this. Thus, for the range < e < 2/3, 
the radial infall model produces flat rotational curves, i.e. it is in accord with the main 
feature of the galactic mass distribution. 

To fit the model to our present galactic halo we must determine appropriate values of 
i?(io) and M(to), where to is the age of the universe. First, we have: 

RitoY 6.52 X lOVar 
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2^2GM(to) h ' 

where h parametrizes the present Hubble rate: Hq = /i 100 km sec~^ Mpc^^. The first 
equality in Eq.(^) follows from the fact that a given shell does not cross any other shell till 
after the moment of its first turnaround and therefore the mass interior to it stays constant 
till then. The second equality in Eq. follows from the assumption Q = 1 which is 
necessary to have a self-similar solution. Second, we match the rotation velocity in the 
model to the observed one for our galaxy: 



where z/(e) is extracted from the numerical solution of the model. Combining Eqs. and 
(I), we have: 

, 1.32 Mpc 

For typical values of e and h, R{to) turns out to be in the 1 to 3 Mpc range; see Table I. This 
is consistent with the value one would infer by taking nearby galaxies as tracers of mass. 



Indeed M31, at a distance of 0.73 Mpc, is falling towards us with a line-of-siglit velocity of 
120 km/s, whereas galaxies at distances exceeding 3 Mpc are receding from us as part of 
the universal Hubble expansion. In fact, turning this around, Eqs. (^) and (|^) could form 
a basis for estimating the age of the universe. 

Figs.l and 2 show the phase-space diagram and the velocity peaks on earth for e = 0.2 
and h = 0.7. Our distance to the galactic center is taken to be rs = 8.5 kpc. The rows 
labeled j = 0.0 in Table I give the density fractions and kinetic energies of the first five 
incoming peaks in the radial infall model. For each incoming peak there is an outgoing peak 
with approximately the same energy and density fraction (see Fig. 2). Note that, because of 
the scale invariance of the model, increasing h at fixed e is equivalent to decreasing R{to) 
keeping the observation radius rs fixed, or to increasing rg at fixed -R(to)- We find that, in 
the radial infall model, the sizes of the two peaks due to particles falling in and out of the 
galaxy for the first time are large, each containing of order 10% of the local halo density for 
e in the standard CDM model inspired range of 0.1 - 0.3. 

However, because the radial infall model neglects the angular momentum that the dark 
matter particles are expected to have, we must question its reliability in this context. In 
particular, since in that model all particles go through the galactic center at each pass, the 
halo density is large there, behaving as phaio(^) ~ l/r^y^ln(l/r). Instead, the actual halo 
mass distribution has the approximate form Phaio('^) = Phaio(0)/(1 + (r/a)^), where a, called 
the core radius, is of order a few kpc 0. Since rg and a have the same order of magnitude, the 
corrections due to angular momentum are not small. Moreover, angular momentum affects 
most the peaks with the highest energies (the ones we are most interested in) because these 
are the ones due to particles which have come from furthest away. 

Fortunately, there is a generalization of the radial infall model which takes angular 
momentum into account while still keeping the model tractable. In this generalization, each 
shell a is divided into subshells labelled by an index k, (k=l, ... , N). The particles in 
a given subshell {a,k) all have the same magnitude lk{<y) of angular momentum. At each 
point on each subshell, the distribution of angular momentum vectors is isotropic about the 
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axis from that point to the galactic center. Thus the spherical symmetry of each subshell is 
maintained in time. Moreover, the evolution is self-similar provided: 

4(a) = ikr^{af/U{a) , (6) 

where the jk are a set of dimensionless numbers characterizing the galaxy's angular mo- 
mentum distribution. In Eq.(^ we are neglecting the small dependence of the turn-around 
radius r*(a) and time t*(a) upon fc, although that neglect is not necessary for self-similarity. 



The model will be described in detail elsewhere [10|. Let us just give here the two equations: 



''''-'^ ^''""m{^\, (7a) 



A<K) = ||:"./°°;i^«(«-;^), Pb) 

and the boundary conditions: Afc(l) = 1, dXkldril) = 0, which determine the model's 



evolution through rk{a,t) = r^{a)Xk{t/t^{a)) and M{r,t) = M{t)M{r / R{t)). In Eq. ([TaD, 
rtk is the mass fraction contributed by the subshell k. In all cases presented here, the jk are 
taken to be distributed according to the density 



exp(-jVjo') • (8) 



We report our results in terms of the average j = -yvrjo/S. Fig.3 shows the rotation curves 
for the case e = 0.2, j = (whose phase-space diagram and velocity peaks are shown in 
Figs.l and 2) and the case e = 0.2, j = 0.2. It shows that the effect of angular momentum 
is to give a core radius to the galactic halo distribution. The definition of core radius we use 
below to report our results is the radius b at which half of the rotation velocity squared is 
due to the halo. For the density profile p(r) = p(0)/(l + (r/a)^), b = 2.33a; of course, our 
density profiles are only qualitatively similar to that one. Fig. 4 shows the velocity peaks for 
the case e = 0.2, j = 0.2, h = 0.7. Table I gives the values of the current turn-around radius 
Rito), the core radius b, the halo density at our location p(rs), and the density fractions 



and energies of the five most energetic incoming peaks for various values of e, j and h. The 
range of e values chosen is motivated by models of large-scale structure formation ^ as well 
as by the flatness of the rotation curves produced. Table I is the summary of our results. 

In conclusion, we have generalized the radial infall model of galactic halo formation to 
include the effect of angular momentum, all the while keeping the model spherically symetric 
and self-similar in its time evolution. The galactic halo distributions we obtain have core 
radii as well as fiat rotation curves. We find that the contribution to the local halo density 
due to particles which are falling in and out of the galaxy for the first time or which have 
passed through the galaxy only a small number of times in the past, and which are therefore 
not thermalized, is rather large, comprising several percent per velocity peak. Finally let us 
emphasize that the peak sizes we obtain are only order of magnitude estimates since they 
are averages over a broad distribution of possible angular momenta. 
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TABLES 

TABLE I. Relative magnitudes Ak and kinetic energies Ek of the first five incoming peaks for 

various vahics of e, j and h. Also shown arc the current turn-around radius R in units of Mpc, the 
core radius b in kpc, and the local density p in 10~^^ g cm~^. The A^. are in percent and the E^ 
are in units of 0.5 x (300 km s~^)^. 
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FIGURES 
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FIG. 1. Phase space of halo dark matter particles at a fixed moment of time with e = 0.2 and 
j = 0. Dotted line corresponds to the Sun's position if h = 0.7. 
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FIG. 2. The spectrum of velocity peaks at the Sun's position for the case e = 0.2, j = and 



h = 0.7. 
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FIG. 3. Rotational curves for the case e = 0.2, with and without angular momentum. 
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FIG. 4. The spectrum of velocity peaks for the case e = 0.2, j = 0.2 and h = 0.7. 
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